Non-Singular Perturbations in a Bouncing Brane Model 
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The question of how perturbations evolve through a bounce in the Cyclic and Ekpyrotic models 
of the Universe remains a topical one. Issues concerning singularities at the background level and 
at the level of perturbation theory continue to be a matter of debate. In this report we hope to 
demonstrate a non singular collision between the boundary branes at the background level, and 
circumstances under which all perturbation variables remain bounded through the collision. 

As expected, we find most collisions to be singular even in the full 5-D formalism, where first 
order perturbation theory breaks down for at least one perturbation variable. Only in the case that 
the boundary branes approach each other with constant velocity shortly before the bounce, can a 
consistent, non singular solution be found. It is then possible to follow the perturbations explicitly 
until the actual collision. In this case, we find that if a scale invariant spectrum developed on the 
hidden brane, it will get transferred to the visible brane during the bounce. 

PACS numbers: 11.25.-w,04.50.+h,98.80.Cq,98.80.-k 



I. INTRODUCTION 

It has recently been proposed that the " Ekpyrotic" sce- 
nario of the Universe 0, or its revised version, the 
Cyclic scenario 3J, could provide an alternative to in- 
flationary cosmology in providing an explanation for the 
observed spectrum of fluctuations which is nearly scale- 
invariant and adiabatic 4]. The effective four space- 
time dimensional picture of the Ekpyrotic scenario cor- 
responds to an initially contracting Universe which un- 
dergoes a bounce and re-emerges in an expanding phase. 
In the four-dimensional picture, the Ekpyrotic scenario 
is closely related to the Prc-Big-Bang scenario 5] in- 
troduced several years earlier j3J|. The Ekpyrotic sce- 
nario, however, is intrinsically five space-time dimen- 
sional. Our Universe is a boundary brane in a five- 
dimensional bulk. The phase of contraction (in the effec- 
tive four-dimensional picture) corresponds to the period 
in which a second brane (a bulk brane in initial scenario 
fil, the second boundary brane in the later proposals 
[2 0) approachcs our brane. The bounce is the time 
of collision, and the phase of expansion corresponds to 
the period when the distance between the branes is once 
again increasing j3(|. The bounce was considered to be 
non-singular from the five-dimensional perspective in the 
initial Ekpyrotic scenario, but singular in the subsequent 
versions. 

It has been shown || that in the contracting phase 
of the Ekpyrotic scenario quantum vacuum fluctuations 
generate a scale-invariant spectrum of adiabatic cos- 
mological fluctuations (more precisely, a scale-invariant 
spectrum of perturbations of the Bardeen variable <&, the 
fluctuation of the metric in longitudinal gauge - see [lj} 
for a comprehensive review of the theory of cosmological 
perturbations, and [ll| for a recent condensed version), 
and that this scale-invariant spectrum would survive the 
bounce (see |l2l for a general criticism of the matching 
conditions of [S| at the bounce). However, if one uses 
Standard matching conditions [Ta PLH ] on a constant field 



hyper-surface to match the fluctuations in the contract- 
ing and expanding phases, then one finds 0] that the 
scale-invariant mode of the contracting phase does not 
couple to the dominant mode of the expanding phase, 
and that the dominant mode in the expanding phase 
will inherit a phenomenologically completely unaccept- 
able n = 3 spectrum (see also [l(| Ht)). Alternatively, 
one finds that even in the contracting phase, the fluctu- 
ation variable C, which is the variable in terms of which 
the fluctuations are well-behaved on super-Hubble scales 
in inflationary cosmology, has a blue n = 3 spectrum |18| . 

The question of how perturbations evolve through the 
bounce in Ekpyrotic and Cyclic cosmologies (in terms 
of the bounce there is no diflerence between the second 
Ekpyrotic scenario of and the cyclic model) is thus 
of central importance. In 0] it has been stressed that 
the spectrum of fluctuations obtained using a matching 
surface analysis depends crucially of the exact choice of 
the matching surface, and in pol | it was pointed out that 
in any coordinate system, one of the pair of canonical 
variables for the fluctuations diverges, rendering the ap- 
plicability of any perturbative analysis questionablc. 

There have been recent attempts to - within an effec- 
tive four-dimensional picture - regularize the bounce and 
follow the fluctuations through the bounce. One attempt 
[2l| was to regularize the bounce using higher derivative 
terms in the gravitational action (similar to what was 
done for the Pre-Big-Bang scenario in with results 

supporting a final n — 3 spectrum for the cosmological 
fluctuations (see, however, the caveats discussed in |23j). 
In contrast, in the work of [2J, |25| it was found that 
when modeling a nonsingular bounce using a spatially 
closed Universe with matter obeying the usual energy 
conditions, the final spectrum can depend crucially on 
the details of the bounce. 

Thcrefore, it seems that only a fully five-dimensional 
analysis of fluctuations can give the complete answer 
as to the final spectrum of perturbations in the Ekpy- 
rotic scenario. However, once again the singularity of 
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the bounce provides problems. A particular matching 
condition, based on the assumption that the bounce in 
the t — y plane (y being the coordinate perpendicular to 
the branes) is a compactified Milne space, was recently 
studied in (2(| (based on earlier work in j2]J - see also 
psj). yielding a scale-invariant spectrum. On the other 
hand, the singularity of the bounce could be more compli- 
cated than a simple Milne space (it could include whiskcr 
regions), and in this case the influence of the whisker re- 
gions on the post-bounce observer region would destroy 
the scale-invariance of the spectrum after the bounce (see 
e.g. HJ) : 

Thus, in this work we will turn to a study of the evo- 
lution of fluctuations in the five-dimcnsional model as- 
suming that the bounce is nonsingular at the level of the 
background cosmology. We will thus assume that the 
collision is a nonsingular event in 5D and only appears 
singular in an effective 4D description. Note, however, 
that we will take the radius of the extra spatial dimen- 
sion to vanish at the time of the bounce, in the spirit of 
the scenarios of 10, but unlikc what occurs in the orig- 
inal scenario of [l|, or in models in which the bounce is 
smoothed out such that the branes never hit. Then one 
can derive, using homogeneity and isotropy in the three 
spatial dimensions, the general form of the metric near 
the bounce [i^J. Under the same assumptions, the field 
content of the bulk can also be computed 0, |3(J . 

Given this background, we will perturb the metric and 
the bulk matter fields. We will then try to find a con- 
sistcnt, nontrivial solution to the perturbed equations of 
motion, such that 

1. The bounce remains nonsingular, that is the Rie- 
mann tensor stays bounded. 

2. Perturbation thcory remains valid, i. e. all per- 
turbed variables stay bounded near the bounce. 

In the majority of cases, it will turn out that there is no 
such solution for the cosmological fluctuations, and thus 
the requircmcnt that the perturbation variables remain 
bounded is untenable. Therefore, the usual techniques of 
evolving perturbations through the bounce are not appli- 
cable. In fact, as a consequence of the divergence of the 
fluctuations, the bounce itself will be a highly inhomoge- 
neous and singular event even when viewed from the 5D 
angle, p ossibly looking like a gas of black holes (see e.g. 
|3ll|33 for recent ideas along these lines). 

The only consistcnt nonsingular background metric 
turns out to be realized if the two branes approach each 
other with constant velocity near the bounce. In that 
case we will be able to explicitly compute how the per- 
turbations evolve during the bounce, so that there are 
no ambiguities regarding matching perturbations before 
and after the bounce. We find in this case that an ini- 
tial bulk scale-invariant spectrum transfers to a scale- 
invariant spectrum on the brane, thus supporting the 
conclusions of 26]. 

The outline of this paper is as follows: After setting up 
a specific realization of a bouncing Universe and review- 



ing the background solution in Section[nJ we will perturb 
the metric and matter fields and identify the scalar de- 
grees of freedom and the resulting equations of motion 
(EOM) f Section 1111(1 . Thereafter we search for a con- 
sistcnt solution, by first idcntifying constraints and then 
working with the equations of motion directly (Section 
II V|) . We conclude with a diseussion of the implications 
of our results for the Ekpyrotic and Cyclic models of the 
Universe (Section IVj l. 

If not stated otherwise the following notational conven- 
tions apply: capital Latin indices run from to 4, small 
ones from 1 to 3 (spatial coordinates), and Greek indices 
from to 3 (time and spatial). Repeated indices are to 
be summed over. A dot denotes a partial derivative with 
respect to the bulk time i, a prime indicates the partial 
derivative with respect to the extra dimension y, and A 
is the 3-d flat space Laplacian. We scale the time t such 
that t = is the time of the collision of the branes. 



II. THE SETUP 

To realize the Cyclic/Ekpyrotic scenario, we will work 
with the aetion of heterotic M-theory [33| 

S = Shet + S BI + S ma tter , (1) 

where Sbi deseribes the brane interaetion that vanishes 
(vanishing at least ~ t near the bounce), S ma tter the 
matter content on the two boundary branes at y = y± = 
and y = yi (taken to be empty before the bounce, and 
containing an ideal fluid or scalar field matter after the 
bounce) and 

S h et = ¥r f d 5 x^g(R - \d A zd A ~ (2) 

2 J M V 2 

^ 1 23 -r- ^ABCDE\ 

— 251 abcdeJ~ j 
2 . 

z=l Jm i 

is a minimal version of the 5D heterotic M-theory aetion, 
with M5 the Planck mass in 5D, R the Ricci scalar in 
5D, S the dilaton, and e a the volume of the Calabi-Yau 
threefold. In addition, Aabcd is the four form gauge 
field with field strength T = dA, eti arc the tensions of 
the visible and hidden branes for i = 1,2 respectively, 
gAB is the metric on M.5 = Mi X Si/Z 2 , hfy is the 
induced metric on the brane, and X^ are the brane 
embedding coordinates. 

The background metric gAB will have positive signa- 
ture, that is 

ds 2 = -n{t,y) 2 dt 2 + a{t,y) 2 S lJ dx i dx 3 + b(t,y) 2 dy 2 .(3) 
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A. Z2-symmetry 

As is already implicit from the form of the action Shet, 
we are assuming a compactification of the five dimen- 
sional theory on the orbifold S\jZ 2 . The visible and the 
hidden brane are located at the orbifold fixed points at 
y = yi and y = y?.-, respectively. It should be noted 
that working on an orbifold yields some additional con- 
straints. For example, the metric has to be Z 2 -even and 
any infinitesimal coordinate transformation £ A must be 
such that £ y (z/ = J/i) = in order to avoid delta-function 
singularities in the metric (see |2J| for a more detailed 
discussion) . 



B. The Einstein equations 

From the metric @ it is straightforward to compute 
the unperturbed Einstein equations, which have to be 
satisfied by the background solution for the metric and 
the fields S and T (see [3(J )■ 



C. Matter equations of motion 

Varying Shet from J2J) with respect to the fields S and 
respectively, yields the EOM 

= U~-^T ABC deT ABCDE 
5! 



+ Y / 5(y-y l )b- 1 6a l e-~, (4) 
= D M (e 2E F MABCD ) 

2 

S%6? ] £ S(y - y^-g)- 1 ^ . (5) 



i=l 



D. Background solution 

The background metric is given by Eq. We write 
each matter field as the sum of a spatially homogeneous 
background field and a field representing small deviations 
around it (and whose spatial average vanishes), that is: 

~(t,x,y) = ^(t,y)+^E(t,x,y) 

ee (°)S + 5S, (6) 

•F ABCDE(t, X, y) = €ABCDE >x,y)) 
= e ABC DE( {0) F + 5F) . (7) 



Note that the Levi-Civita tensor implicitly contains the 
appropriate metric factor so that it is a bona-fide tensor: 

CABCDE '■— \f—g^ABCDE- 

One can derive the general form of the unperturbed 
metric and the unperturbed matter fields given above. 



Assuming non-singularity, then near the bounce the Solu- 
tions can be written in terms of a power series expansion 
in t (see @,|23)- O ne then takes the following steps: 

1. One demands that the bounce is non singular in 
5D, i. e. that the Riemann tensor ^Rabcd, the 
Einstein tensor (°'G A B and the energy momentum 
tensor remain bounded near the bounce. 

2. One then solves the EOM of the matter fields as 
well as the Einstein equations for the unperturbed 
metric for small t. 

3. One enforces the matching conditions arising from 
delta functions, i. e. one has to match the bulk solu- 
tion to the boundary values on the branes at y = y\ 
and y = y 2 . 

This finally yields the following possible forms for the 
metric: 



b(t, y) = Bt« + Bb k+1 (y)t k+1 + 0{t^ z ) 



ik+2\ 



(8) 



2fc-3 

n(t, y) = 1 + Niyf 

i— k 



+n 2k - 1 (y)t 2k - 1 + 0(?«), 



,2k\ 



(9) 



k — l 



2k-2 



a(t, y) = 1 + V A t t' + V (AiBy + A,)f 



i=0 



B 2 A X 



y 2 + A 2 k-iy + A 



2k-l 



2fe-l 



(10) 



2 

+0(t 2k ), 
with k > 3 and 

n' 2 , k _ 1 (y)-B 2 k(k + l)b k+1 (y) = 0, (11) 

and Ai, A4, B and Ni are constants. The Aj may be zero 
and N k = A k = ~a 1 e~~°/2. 

The other allowed metric with k = 1 is [ijj 

b(t,y) = Bt + Bb 2 (y)t 2 + 0(i 3 ), (12) 
n(t,y) = l + m(y)t + 0{t 2 ), (13) 

a(t,y) = 1+ (A 1 smh(By)+A 1 coah{By)jt + 0{t 2 ) 
= l + ai (y)t + 0(t 2 ), (14) 
where n\ and b 2 are related via 



n x (y) n'{{y) 



(15) 



and the boundary conditions are n^(t/j) = a'\(yi) — 
— Ba\e~"° j2. A mistake in |3(j regarding the case k = 1 
was corrected in |35j . Before the bounce k = 1 and k > 3 
are possible, after the bounce only the k = 1 case is al- 
lowed. B should have a different value before and after 
the bounce. 
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Note that demanding in addition to the above bound- 
edness requirements a "no flow condition" (no energy is 
allowed to flow off the brane), then the single allowed 
metric in the case k = 1 after the bounce H12|) - (|14fl i s 
also ruled out |2jj. However, since in the Ekpyrotic and 
Cyclic scenarios one assumes that energy can flow from 
the bulk to the brane, it is not reasonable to demand this 
additional condition. 

For the field one gets by using its equation of mo- 
tion and the junction conditions 

fc-i 

= ^H^ + (3a ie - Ho 5j/ + H fc )t fe 

t=0 

+0(t k+1 ) (16) 

for k > 3 and 

{0} E-{t,y) = E + (Xcosh(By)+Xsmh(Byfjt 

+0{t 2 ) (17) 

for k = 1, where A, A and Hj are constants. The field 
is then given by 

(°^,t)- ai e- 2<0)H ^. (18) 

III. SCALAR PERTURBATIONS 

We will now perturb the metric focusing only on scalar 
metric perturbations. The most general scalar perturba- 
tion in "generalized" longitudinal gauge can be character- 
ized by four scalar functions, <I>, 'f, W and T [34]]. These 
functions can be viewed as a basis of gauge invariant vari- 
ables. In this gauge, the metric is given by 

ds 2 = -n{t,y) 2 (l-2<l>(t,x,y))dt 2 

+a(t,y) 2 (l - 2^{t 1 x 1 y))8 i: jdx l dx j 
+b(t,y) 2 (l + 2T(t,x,y))dy 2 
-2n{t,y) 2 W(t,x,y)dtdy, (19) 

where the signs in front of the perturbations are a mere 
convention. 

Two more scalar degrees of freedom (matter degrees of 
freedom) were already introduced in © and J7J). They 
are <5S(i, x, y) and S J- (t, x, y). 

A. The brane embedding and longitudinal gauge 

It should be noted that we do not perturb the em- 
bedding of the boundary branes - they are still localized 
at y = yi = Q and y = y 2 . We are entitled to do this 
without loss of generality, since one can induce local fluc- 
tuations in the radius by allowing 355 to fluctuate, rathcr 
than by making use of fluctuations in the coordinate in- 
terval. However there are furthcr compelling reasons for 



this approach: perturbing the embedding would render 
the collision into a series of local events in time, thus 
complicating the situation a lot. This approach of not 
considcring perturbations in the position of the branes 
is common in the literature (see e. g. |2(|). However, 
one should note that due to the possible divergence of 
the metric at the bounce, in general a large coordinate 
transformation would be necessary to transform small 
perturbations in the coordinate position of the brane 
away. This could induce further changes in the spectrum 
of perturbations (keep in mind, that the gauge invariant 
variables introduced in H19|) are only gauge invariant un- 
der infinitesimal coordinate transformations). Since we 
are only interested in the physics of how a given per- 
turbation evolves through the global bounce (and not in 
exploring new ways for how the spectrum could be gen- 
erated at an inhomogeneous bounce) our simplified setup 
is satisfaetory. 

Note also that the dynamics of the collision is fully 
includcd in the metric - the brane coordinates are 
time-independent but the physical distance between the 
branes does go to zero when t approaches zero. 

One may wonder if keeping the branes at fbced posi- 
tions is consistent with working in the generalized longi- 
tudinal gauge. First note that the branes sit at the orb- 
ifold fbced points - demanding that the metric may not 
include delta funetion singularities requires an infinites- 
imal coordinate transformation £ to be continuous ev- 
erywhere, specifically £ y (y — y t ) = has to hold at the 
brane positions j34|. This ensures that the coordinate 
transformation necessary to go to the longitudinal gauge 
does not perturb the brane position. 

It is elear that this situation depends on our choice to 
work with branes localized at orbifold fixcd planes. As 
diseussed e. g. in [3(|, in a more general brane world set- 
ting one needs to include an extra degree of freedom for 
scalar metric fluctuations. In our setting, going to longi- 
tudinal gauge and keeping the brane positions fbced does 
not introduce any new restrietions on the perturbations 
considered. 



B. The perturbed Einstein equations 

To write down the Einstein equations we need both 
the perturbed energy-momentum tensor ^'T^ and the 
perturbed Einstein tensor ^G A B , computed up to first 
order in the perturbation variables. The latter is straight- 
forward to compute from l|19|) and reads in component 
form (no summation convention where repeated indices 
appear) 

= -l^.^ + s-r] , (20) 



■5 



(i) G t 



(l)QXi 



(l) G x 




ab a 



n 2 <9„ W 



4 [d xt d Xi - A] (-r + * + + 



a' b' n' , , 

3 — H \-Oy\dy 

a b n ' ■ 



2 I b dd d& a 2 nb 'a b . d „ \ 1 f b' n' a' „ \ » 
— 7-2— + 2- + — - — + 2- + -9 4 + -9 t +72 — r + 2- + 2-+a„ U, 
n 2 \ o na ab a z nb a 2b a I b z \ b n a J 



— ( 2 2 

6 2 V a ab 



n'b' 
nb 



n" n' a' 
— + 2 

n na 



a' 



(L 



—i 

2n 



„a n „b _ \ _ 

2 + 2 7 + aU 



2 
"b 2 



n' „ d' fta' 
— + 2- + 

n a na 



n6' dn' 
2nb an 



bn' b'a ba' 
2bn ba ba 



2b 



-U 

n 



2^ 



9t 
2 



(21) 



o 2 



a' 
a 



n' 
n 



n 



a 



+ 



1 a 6 
cr n A 



b 2 
d 



er \ a n 



2 o T 1 / a' n' n \ ^ 1 a' n' \ ^ 1 b a h „ 



2fe 2 



2 . 3 /_ n , d\ _ 3 / „a' n' \ _ 
-^A + — 3 t - - + 4- 0t + - -2 A. 



a 



n 



b 2 



b 

b a 

/2 



W 



T 6 ( a 



a z 



an 



r 



d nd a. d \ 3 a' 

— + _ + ^T^ ~ 12 — 9 2/ 

a z na a 2a / tr a 



$ - 



6 2 



na' „a' d a'6 d' a' \ 
+ 2— - — + - + -9 t . 

na a z ab a a 



(22) 



(23) 
(24) 



(25) 
(26) 



The right hand sidc of thc Einstcin cquations is given by 
the perturbed energy-momcntum tensor, which can be 
derived in a straightforward way by varying the action 



with respect to the metric and keeping only first order 
contributions. One finds that 



mi 



% c (0 *d c 5z + ld A s~d B + \d A ^hd B SE + s A l±-e 2 <0>s ^t CD efg ^t cdefg 



5 1=1 



1 sp s (y-yi) 

b( yi ) 



^94 , . . 



(si 



3 ai Ml-T A B 



2 5! 

S(y-yi) x 
-btmT 5gBfiT ^ 



6g Bf3 T A f + 6 A e- H 3a 4 M 5 3 <5S + 5B/3 oT ( 



(27) 



where the last sum is only present if A ^ 4 and B^4. ^From the a^-a^-component of the perturbed Einstcin 

equations it follows that the the a^-a^-component of 
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^G A B vanishes, as can be seen from (|27|) . Thus from 
(l2T)jl we have 



[* + $ - r] = o . 



(28) 



By exploiting rotational invariance in the three spatial 
directions, and requiring that the perturbations remain 
bounded everywhere and have a vanishing spatial average 
one can deduce that Y = ^ + <f> up to a function of t and 
y Hl|. Thus, 

T{t, x, y) = $(t, x, y) + *(t, f, y) + /(i, y) , (29) 

where f {y, t) is a small, so far undetermined function 
whose average over the extra dimension has to vanish. 
Thus, fbr any spatial Fourier mode except the zero mode 
we may write 



T(t,q,y) = $(t,q,y) + V(t,q,y), 



(30) 



and the perturbed metric l|19|) . One gets 



n z 
1 




9,1 + - - -r + 3- 5, 



(1) D = 2n 2 $<9 2 - 2a 2 *A + 26 2 r<9 2 - 2n 2 Wd y d t 




h a b \ 
-2- + 6- + 2- U> + $-3* 
n a o / 



6' 



FF + VF' 



(34) 



2— - 2^- + 6— ) r + r' + $' + 3*' 



bn a\ 
-r + - + 3- FF + VF 

bn a 



1 3 

_^($ ii + * ii -r i )5 Xi) 



(35) 



but the zero mode still contains the function f (t, y). 
Since we are doing linear perturbation theory, all Fourier 
modes decouple and thus f (t, y) is of no importance 
for the spectrum of metric fluctuations. However, it 
may be important for the evolution on the brane, since 
f {t, y = j/i) may look like a small, time-varying cosmolog- 
ical constant on the brane. This subtlety is often glanced 
over in the literature 1341. 



where A is the usual spatial Laplacian. 

Note that l|32|) is in fact a algebraic equation relating 
S J- to (5S, since a simple integrations yields 



8 T = -2 (°JF5S , 



(36) 



where the integration constant was set to zero since this 
is the conclusion one would arrive at from solving for the 
gauge field in terms of the dilaton through the equations 
of motion 301. 



C. The perturbed matter equations of motion 

To get the perturbed equations of motion we have to 
plug in the perturbed metric from i|19[l into 10} and JSJl, 
use the definitions © and J7J), and expand the resulting 
equations up to first order in the perturbations. For the 
smooth componcnt of the fields we obtain 



+4 (SE^T 2 + (°^)e 2(0)H 
= dJ\^T25E + 5T\ e 2< ° )H 



(31) 
(32) 



We will return to the delta function part of (|31|l in seetion 
IIVA3I where we will extract constraints from it. (°tl and 
( x tl are straightforward to compute using the defmition 



□ 



-AB- 



9 BC T A B C dA 



g^d A d B - g»"Ti 
(o)D + (i)n + 0(perturbations 2 ) (33) 



IV. IN SEARCH OF A CONSISTENT 
SOLUTION 

Now that we have introduced our notation and 
equipped ourselves with all the relevant equations of mo- 
tion we will try to find a self consistent, nontrivial solu- 
tion. Let us first make a general ansatz for all perturbed 
quantities and try to exploit all constraints before plug- 
ging the ansatz into the respective equations. 

There are three main sources of constraints: 

1. In order to be consistent with our background so- 
lution, we must demand that the perturbed Rie- 
mann tensor ^Rabcd remains bounded near the 
bounce. 

2. The perturbed energy-momentum tensor '''Tj 
should remain bounded. 

3. The delta function parts of the EOM have to be 
matehed with possible divergent parts originating 
from derivatives with respect to the extra dimen- 
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Assuming that the fields describing the fluctuations do 
not diverge at the bounce t = 0, we write down a general 
power law expansion for them (following the philosophy 
that was used for the background fields). We will then 
check as to whether nontrivial solutions are possible. If 
not, we can then conclude that the perturbations diverge 
at the bounce. Thus, we make the ansatz 



$(t,x, y) 

y) 

W{t,x, y) 
SE(t,x, y) 



<p (x,y) + <p s (x,y)t s +0(t s+1 ), (37) 

Mx, y) + M& v)) tq + o( tq+1 ) . (38) 

w (x,y) + w l (x,y)t l +0(t l+1 ), (39) 
y) + y)t m + 0(t m+1 ) , (40) 



where s,q,l and m are the positive exponents of the lead- 
ing order contributions in t. 

If we want the bulk perturbations to be able to seed 
nontrivial branc fluctuations during the bounce, we re- 
quire the additional condition ipo(x,y) ^ and/or 
il>o(x,y) if k > 3. The possible non-trivial case that 
the perturbations vanish but their time derivative does 
not will be ruled out shortly, since it turns out that q > k 
and s > k, with k > 3. 

Note that if this ansatz is consistent and nontrivial 
near the bounce, any perturbation variable evaluated on 
the visible brane at the bounce is trivially determined by 
its bulk value shortly before the bounce. 

This is consistent with the physical picture of the 
Ekpyrotic/Cyclic scenarios: Perturbations are generated 
as quantum vacuum fluctuations of the bulk scalar modes 
well before the collision of the branes (pictorially, they 
can be represented as fluctuations in the distance be- 
tween the branes). They evolve as the two branes ap- 
proach each other until shortly before the bounce. The 
key question is how the initial spectrum of fluctuations 
in the "contracting" phase (the phase during which the 
branes are approaching each other) is transferred to the 
spectrum of perturbations on the brane immediately af- 
ter the bounce, assuming that there are no pre-existing 
fluctuations on the brane. In the Ekpyrotic/Cyclic sce- 
narios, the background evolution is chosen such that this 
bulk spectrum is scale-invariant 42]. Our analysis, how- 
ever, is independent of the origin and shape of the initial 
bulk spectrum. 



Before the bounce 



There are two possibilities k — 1 and k > 3. We 
will first work out some general constraints before we 
study these two cases separately. In addition we will 
only focus on the spatially inhomogeneous modes of the 
perturbations, so all results should be taken modulo a 
possible non-vanishing zero Fourier mode. 



1. Boundedness o} the Riemann tensor 

Demanding that potentially divergent terms in 
^Rabcd get canceled, we obtain the following con- 
straints on the exponcnts: 



/ > fc + 1, 
s > k , 
q > k. 



(41) 
(42) 
(43) 



Furthermore, the functions in our ansatz l|37(l - (|4()|l are 
constrained to 



wo{x,y) = 0. 



(44) 
(45) 
(46) 



Deriving this is a straightforward exercise - one simply 
writes down all components of the Riemann tensor to all 
orders in t until C(t _1 ) and demands that they vanish. 

To be specific, let us retrace the main steps leading to 
the constraints in the case that k > 3: ^Rtx 3 x 3 y yields 
to order O^' 1 ) 



% + —w = . 



(47) 



On the other hand, ^■ 1 'Rt XlXl x i yields to order 0(t ■ 7 ) 
AA k 



B 2 



-d Xi w = . 



(48) 



so that wq can not depend on x. Therefore wq only cou- 
ples to the zero Fourier modes, which we are not inter- 
ested in. Thus, we may set wq = 0, and from (|47|l we get 
consequently ipo(x,y) = ipo{x). Stepping up in orders of 
t in ^Rtx 1 x 1 x i to (9(t -J ), where I < j < k, we obtain 



d Xi Wj = , 



(49) 



so that l > k. Now, looking at ^Rt Xl tx! to order <D{t~ k ) 
we get, using the previous results, ipo(x,y) = (po{x). 
Turning to ^R Xl x 2 x 1 x 2 we obtain to leading order 



2A A k ^ - 2A 1 A k w l t l = 0(t k ) , 



(50) 



so that q > k has to hold. As a consequence, it follows 
from ( 1) R t x 1 tx 1 



<p' s t s = 0{t k ) . 



(51) 



so that s > k has to hold too. Finally, if we consider the 
2/-Xi-Einstein equation to order 0(t~( k+lS> ), we see that 



d Xi w k = 0(t k+1 ) 



(52) 



must be satisfied so that l > k + 1 - all in all we have 
derived the claimed constraints in the case k > 3. No 
other component of the Riemann tensor gives additional 
constraints. The arguments for k = 1 are similar. 
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2. Boundedness of the energy momentum tensor 
Looking at 

1 



we see immcdiately that 



m > k 



(53) 



(54) 
(55) 



has to hold in order for the energy-momentum tensor to 
remain bounded. The other components do not give any 
additional constraints. 



3. Matching to the branes 

The next conditions we should look at originate from 
matching <5-function parts in the Einstein equations and 
the EOM of the field SE. Note that the background So- 
lutions already satisfy similar constraints including only 
unperturbed quantities. 

Let us first look at the Einstein tensor: the only pos- 
sible divergent terms consist of second order derivatives 
with respect to y. Thus, from l(2T j) -ip7 )l we see that only 
the t-t and X\-Xi components of the Einstein tensor con- 
tain such terms. These parts are, specifically, 



t\s 



(l) ff 



b 2 b 2 a 1 

2 

~52 





(la" 


n"\ 1 




h— r 




V a ' 





(56) 
.(57) 



Note that the summation convention is not used in the 
above expression for the x%-Xi component. Let us intro- 
duce some notation |37| : Splitting X" into a continuous 
part X" and a delta function part multiplied by the jump 
in the first derivative [X'], that is, e. g. 



$" = d" + [$'] 6(y - Vi 



(58) 



Now, exploiting the Z2-symmetry of the bulk we know 
the jump of the perturbation variables. For example, in 
the case of the brane at the origin we know 



<Z>'(t,x,y = 0^ 
*'(*,£, y = 0^ 
T'(t,x,y = 0^ 
W'(t,x,y= O 4 



= -$%2,y = 0-), 

= -*'{t,x,y = 0-), 

= -T'{t,x,y = 0-), 

= W'(t,x,y = Q-). 



Thus, e. g. for $' at y = we have 
[$'] 



<P'(t,x,y = + )-<S>'(t,x,y = 0-) 
23>'(t,x,y = 0+). 



(59) 
(60) 
(61) 
(62) 



(63) 



Note that a, b and n behave like $ at the brane positions 
under the Z 2 -symmetry. With the perturbed energy- 
momentum tensor of the form (j27(l with Q56|l and Ij57(l 
we get our first two constraints 



6_M 

6 2 a 



4( 2 M + E 

b \ a n 



r. 



which have to be satisfied on the branes at y — and 
y = J/2- We have made use of the assumption that there 
is no matter present on the branes before the collision. 
The physical picture is, of course, that matter is created 
during the collision by the transfer of kinetic energy of the 
hidden brane to the visible brane, but for our purposes it 
will be satisfactory to simply assume a sudden creation of 
matter at t — 0, so that we have some ideal fluid present 
on the visible brane for t > (naturally, this assumption 
does effect the matching conditions for positive t). 

A more convenient form of these two constraints is ob- 
tained by taking the difference and the sum, yielding 



O(t) 



[$'] - [*'] - 2 



a 



n 



(64) 



and 
6 



i (T + 5E) + 0(t) = -- (65) 




+2 5 



We obtain one more constraint from the delta function 
part of the equation of motion for SE in Q, expanded 
to first order in the perturbations. One gets, e. g., for the 
brane at y = 



= o(t). 



(66) 



4- The Einstein Equations for k > 3 

In case q > 3, it follows from the junction conditions 
(UU and that 



s > 3, 
(fo(x) = -VoO?) ■ 



(67) 
(68) 



In case when q = 3, all we can deduce from the junction 
conditions alone is 

s = 3, (69) 
^'k{x,y = Vi) = ^'s{x,y = yi) 

~ (yoffl + MZ))- ( 7 °) 
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Now, looking at the Xi-t Einstein equation (|22fl to lead- 
ing order in t, we see that w' k+l (x,y) = 0. Using this 
knowledge in the y-x% equation (J2U, we get 



w k +i(x,y) ~ {<po(x) + tpo(x)) . 



(71) 



Working out the y-t equation (|23|1 to order O {t k ) and 
evaluating it at the brane we obtain, after using (|70|l and 

CD. 

ipo{x) = -ipa{x). (72) 
Thus, in any case we can deduce l > k + 1 as well as 

ip (x) = -ipo(x) (73) 

and 

ifg=fc=> f q (x,y = yi ) = 0, (74) 

ifs=fc=> <p' a (x,y = yi) =0. (75) 

Going back to the junction condition for <5S, it follows 
that 



if m 



Cn&V = Vi) =°- 



(76) 



Equippcd with this knowledge, it can be seen that the 
leading order contribution to the y-y equation (|27|l is of 
order O(t ) and takes the form 



Aip (x) + Cip a (x) = 0, 



(77) 



where C is a constant given by —6A 2 — I2A2. Thus, 
only one specific Fourier mode of ipo could be consistently 
transfcrred to the visible brane. This certainly is incon- 
sistcnt with obtaining a scale invariant spectrum on the 
brane. Thus, we will not further consider this possibility. 
Therefore, we conclude that only the trivial solution 



ip Q (x) = 0, 
M*) = 



(78) 
(79) 



is consistent with our model if k > 3. This is our first 
main result. Its implications will be discussed in section 

m 



5. The Einstein Equations for k = 1 

Let us consider the case k = 1, that is the metric 1121) - 
ljTl)l. The Xi-t equation evaluated to order 0(t~ r ), 
yields 



1po(x) = -2<y3o(^) • 



(80) 



Since the junction conditions l|64|) and i|66|) imply 
ipo(x) = —ipo(x) if s > 1 and/or q > 1, we must have 
s = q = 1 and the boundary conditions 

Vi&y = Vi) = a'i(y = Vi) ( Po(x), (81) 
$'i{x,y = Vi) = a'iiv = Vi)'fio{x), (82) 
£i(x,y = Vi) = 6ai(y = y t )ip (x) ■ (83) 



The y-Xi equation (|25|l to order O (t 1 ) yields 



2ip[ + ip[ - 3a[ip Q + -w 2 = . 



(84) 



Evaluating this at the brane gives us the boundary con- 
dition 



w 2 (x,y = j/i) = 0. 



(85) 



Anothcr equation of motion follows from the t-t equation 
(ED to order O^ 1 ) 



B 2 



1>i 



which is solved by 

i>i($>y) = ai(y)<po(x) 



(86) 



(87) 



where we made use of (|82|) . The Xi~t equation l|22ll to 
order O(t ) gives yet another equation of motion 



ipi +a 1 ip Q + -^ w 2 = 0, 
which can be combined with l|84|) to yield 



B 2 



9(fi + lOai^o , 



w< 2 r 

-^2 = 9w 2 - 120^0 • 

These are solved by 

5 3 ~ 

fi(x,y) = -~ai(y)<p @) + -Ai<p (x) 



(89) 
(90) 

(91) 



, ' 1 '. COSh(3 ^ 2) cosh(3i? y ) + S mh(3By) ) , 



smh(3By 2 ) 



w 2 (x,y) = 2 BAl 



+ 



cosh(3By) — cosh(By) 
1 — cosh(3_B?/2) 



(92) 



sinh(3By2) 
1 - cosh(y 2 -B) 
sinh(j/2-B) 



■ sinh(3By) 



smh(yB) 



<Po(x), 



where we made use of the boundary conditions l|81|l and 
(|85|l . For £1 we get from (I31|) . to leading order in t, 



— = £ x + I2anp , 



(93) 



with the solution 

ti{x,y) - 



3(A 1 + A 1 )[By + C\--)e B y (94) 



3(A 1 +A 1 )(By + C 2 + 



-By 



Vo{x) , 
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where C\ and C 2 are constants chosen such that l|83|) is 
satisfied. A nontrivial consistency check for these Solu- 
tions is, e. g., the Xi-Xi equation l|24|l evaluated to order 
©(t" 1 ), which is indeed satisfied. 

Now we have to look at all other Einstein equations 
up to order OliP), all of which include higher order per- 
turbations. In order for a scale-invariant spectrum to be 
consistent, the terms ~ Ay>o in l)27|l and 1|21[) should get 
canceled for all y. We can achieve this if we write 

ip2(x,y) = fa(x)ip 2 (y) + ^2(x) (95) 

and demand that 

Mx) = ^Aip (x). (96) 

Combining to order O^ 1 ) with (|2TJ and J23 to or- 
der O(t ) gives three equations of motion for 02(2/), ni(y) 
and %p 2 {y). Together with ((To^l these could be solved. 
Since a 2 is a solution to the background equations involv- 
ing the brane interaction, one could in principle extract 
information about the brane interaction by demanding 
consistency in the perturbed metric beyond leading or- 
der. Since we did not specify the brane interaction so far, 
we shall be satisfied with consistency to leading order. 

The remaining three Einstein equations i|23|) , (|25|1 and 
(|24|l . all evaluated to order O(t ), involve six new degrees 
of freedom: ^3, ip 2 , W3, n 2 , 63 and 03. These define tp3, ip 2 
and u>3 in terms of higher order terms of the metric, all 
of which require a better understanding of the brane in- 
teraction. 

All that is left to do now is to check if we can impose 
the initial conditions 

<Po(£) = ~<Pi(x, V = 0)*, , (97) 

M g ) = -M^y = o)u, (98) 

such that at t — ti perturbations <& and 'f 011 the visible 
brane vanish up to order O (t?) (they do not vanish on 
the hiddcn brane). We can arrange this if we assume 
By 2 » 1, -aie~ So > 8 - then U > -1/2 and 

®{x,y = y 2 ,t = U) = ^ (S)+O{tj), (99) 
V(x, y = y 2 ,t = U) = -4p {x) + 0{t}) . (100) 

It should be noted that in order for this to be possible 
we need the visible brane to be a negative tension brane. 
Since at t = the perturbations on the visible brane are 
given by $ — ipo and 4" = — 2(fQ, the above two equa- 
tions tell us how to match perturbations shortly before 
the bounce in the bulk (and on the hidden brane) to 
perturbations on the visible brane at the bounce. 

The only change to our setup after the bounce is a 
change in the junction conditions in section IIV A 31 due 
to the matter which has been generated by the bounce on 
the visible brane. Therefore, a consistent solution after 
the bounce could be found along similar lines as above. 



V. CONSEQUENCES FOR THE EKPYROTIC 
AND CYCLIC SCENARIOS 

We have seen that for k > 3 the only solution to the 
perturbation equations that is bounded near the bounce 
and does not yield a singularity in the Riemann tensor 
is the trivial solution l|78|) and l|79|) - the solution that 
has vanishing metric perturbations on the branes at the 
bounce. 

Since there have to be some perturbations on the vis- 
ible brane after the bounce if the Ekpyrotic Universe is 
to describe our Universe, we are led to conclude at least 
one of the following 

1. The background Riemann tensor and/or energy- 
momentum tensor does not remain bounded at the 
bounce and thus we have a singularity even from 
the 5D point of view. 

2. At least one perturbation variable exhibits a grow- 
ing mode ~ t~ c where c > 0, and thus perturbation 
theory breaks down as we approach the bounce. 

In either case the bounce can not be examined using Stan- 
dard methods. One has to specify the type of the sin- 
gularity and find a way to match perturbations shortly 
before and shortly after the bounce. Unfortunately the 
breakdown of first order perturbation theory renders the 
whole problem highly gauge dependent. The most recent 
attempt addressing this problem j2(J assumes a Milne- 
like Universe and works in a specific gauge, but the ro- 
bustness of the results towards changes in the specific 
assumptions concerning the background and the match- 
ing conditions is subject of an ongoing debate. 

However, in the case that the two boundary branes ap- 
proach cach other with constant velocity near the bounce, 
that is k = 1, a consistent solution of the perturbed Ein- 
stein equations can be found, as was done explicitly in 
section HY A 51 In this case, one can easily follow the per- 
turbations until the actual collision. Furthermore, if a 
scale invariant spectrum in the bulk was generated while 
the hidden brane was approaching the visible brane, then 
it will indeed be transferred to the visible brane during 
the collision. Since the bounce is non-singular and all 
perturbations remain small, one needs neither to worry 
about higher curvature correetions nor gauge ambigui- 
ties. In order to see if this single case is actually realized 
one needs a better understanding of the brane interac- 
tion and check consistency beyond leading order in the 
perturbed metric. 

The limitations of our conclusions lie in the specific 
realization of the model used. Our choice was guided by 
the principle of simplicity, apparent e. g. in the choice 
of the bulk fields. Furthermore we did not specify the 
brane interaction - we only assumed that it vanishes at 
least like t during the bounce. We showed that in this 
simplified setup a consistent realization of the Ekpyrotic 
or Cyclic seenarios including perturbations is possible. 



11 



Acknowledgments 

We wish to thank Syksy Rasanen for useful correspon- 
dence, and Paul Steinhardt for comments on the first 
version of this paper. One of us (SP) wishes to acknowl- 
edge Jean-Luc Lehners at Imperial College for discussion 



concerning brane perturbations, and the Imperial College 
Theory Group for the use of their facilities in the course 
of the work on this paper. RB wishes to acknowledge 
the hospitality of the McGill Physics Department. This 
work was supported in part by the U.S. Department of 
Energy under Contract DE-FG02-91ER40688, TASK A. 



J. Khoury, B. A. Ovrat, P. J. Steinhardt and N. Turok, 
Phys. Rev. D 64, 123522 (2001). 

J. Khoury, B. A. Ovrut, N. Seiberg, P. J. Steinhardt and 

N. Turok, Phys. Rev. D 65, 086007 (2002). 

P. J. Steinhardt and N. Turok, Phys. Rev. D 65, 126003 

(2002). 

D. N. Spergel et al, Astrophys. J. Suppl. 148, 175 (2003). 
M. Gasperini and G. Veneziano, Astropart. Phys. 1, 317 
(1993). 

R. Kallosh, L. Kofman and A. Linde, Phys. Rev. D 64, 
123523 (2001). 

R. Kallosh, L. Kofman, A. Linde and A. Tseytlin, Phys. 
Rev. D 64, 123524 (2001). 

K. Enqvist, E. Keski-Vakkuri and S. Rasanen, Nucl. 
Phys. B 614, 388 (2001). 

J. Khoury, B. A. Ovrut, P. J. Steinhardt and N. Turok, 
Phys. Rev. D 66, 046005 (2002). 

V. F. Mukhanov, H. A. Feldman and R. H. Branden- 

berger, Phys. Rept. 215, 203 (1992). 

R. H. Brandenberger, hep-th/0306071 

P. Peter, .1. Martili. N. Pinto-Neto and D. .1. Schwarz, 

hep-th/0204227 

J. C. Hwang and E. T. Vishniac, Astrophys. J. 382, 363 
(1991). 

N. Deruelle and V. F. Mukhanov, Phys. Rev. D 52, 5549 
(1995). 

R. Brandenberger and F. Finelli, JHEP 0111, 056 

(2001) . 

J. C. Hwang, Phys. Rev. D 65, 063514 (2002). 

S. Tsujikawa, Phys. Lett. B 526, 179 (2002). 

D. H. Lyth, Phys. Lett. B 524, 1 (2002). 

R. Durrer and F. Vernizzi, Phys. Rev. D 66, 083503 

(2002) . 

D. H. Lyth, Phys. Lett. B 526, 173 (2002). 

S. Tsujikawa, R. Brandenberger and F. Finelli, Phys. 

Rev. D 66, 083513 (2002). 

C. Cartier, J. C. Hwang and E. J. Copeland, Phys. Rev. 
D 64, 103504 (2001). 

C. Cartier, R. Durrer and E. J. Copeland, Phys. Rev. D 
67, 103517 (2003). 

P. Peter and N. Pinto-Neto, Phys. Rev. D 66, 063509 
(2002). 

J. Martin and P. Peter, Phys. Rev. D 68, 103517 (2003). 



[26 

[27 

[28 

[29 

[30 
[31 

[32 
[33 

[34 

[35 
[36 

[37 

[38 



[39] 
[40] 

[41] 



[42] 



A. J. Tolley, N. Turok and P. J. Steinhardt, 
hep-th/0306109 

A. J. Tolley and N. Turok, Phys. Rev. D 66, 106005 
(2002). 

L. Cornalba and M. S. Costa, Fortsch. Phys. 52, 145-199 
(2004). 

B. Craps and B. A. Ovrut, Phys. Rev. D 69, 066001 
(2004). 

S. Rasanen, Nucl. Phys. B 626, 183 (2002). 

T. Banks and W. Fischler, hep-th/0111142 

T. Banks and W. Fischler, |hep-th/0212113| 

G. Veneziano, JCAP 0403, 004 (2004). 

P. Horava and E. Witten, Nucl. Phys. B 460, 506 (1996); 

P. Horava and E. Witten, Nucl. Phys. B 475, 94 (1996). 

C. van de Bruck, M. Dorca, R. H. Brandenberger and 
A. Lukas, P hys. Rev. D 6 2, 123515 (2000). 

S. Rasanen, astro-ph/0208282 

M. Dorca and C. van de Bruck, Nucl. Phys. B 605, 215 
(2001). 

P. Binetruy, C. Deffayet and D. Langlois, Nucl. Phys. B 
565, 269 (2000). 

There are, however differences both in terms of the parti- 
cle physics motivation for the scenario - in the Ekpyrotic 
scenario the string coupling goes to zero at the bounce, 
whereas in the Pre-Big-Bang scenario the bounce occurs 
in the regime of strong string coupling - and also in the 
details of the dynamics - the rate of contraction is differ- 
ent. 

See 0, [7, 8] for some criticism of the basic scenario. 
Note that this metric does not correspond to a Milne like 
bounce, since n\ and ai depend on y. 
One can write the perturbed metric in a rotated coordi- 
nate frame. The infinitesimal perturbation variables do 
not change, and hence we can conclude that all second 
order 4D spatial derivatives of the off-diagonal Einstein 
tensor vanish. 

In these scenarios, the mode of $ with a scale-invariant 
spectrum is growing in the contracting phase so that at 
the time of the bounce the magnitude of the fluctuations 
is much larger than the magnitude of any pre-existing 
quantum vacuum fluctuations on the brane. 



